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Abstract 



A method of reducing general quaternion functions of first degree, i.e., 
linear quaternion functions, to quaternary canonical form is given. Linear 
quaternion functions, once reduced to canonical form, can be maintained 
in this form under functional composition. Furthermore, the composi- 
tion operation is symbolically identical to quaternion multiplication, mak- 
ing manipulation and reduction of systems of linear quaternion functions 
straight forward. 

1 Introduction 

Real linear functions take the monomial form: f(x) — mi, where x, m G K. 
By contrast the general quaternion linear function has a multi-nomial form: 



where all factors are quaternion valued, i.e., q, m p , n p € H. The number of terms 
in the summation can be arbitrarily large due to quaternion multiplication being 
non-commutative. 

Unlike Hardy [TJ p. 199-200] who studied the inversion of these functions 
or Hathaway [5[ chp. 4] who was interested in their properties, this work ad- 
dresses the reduction of systems or networks of such functions, formed by scalar- 
weighted sums and functional composition. The number of monomial terms in 
such systems quickly become intractable. 

Herein is a method for reducing such systems in an efficient and transparent 
manner. Two reduction methods are key to achieving this goal. First is conver- 
sion of general linear quaternion functions, as given in equation |T]), to at most 
a four-tuple form. It will be shown that equation ([1]) can always be reduced to 
the form: 



p 




(1) 



f{q)=Aq + Bqi + Cqj + Dqk, 



(2) 



where A, B,C, D G H and i,j, k are the standard hypercomplex operators. Lin- 
ear quaternion functions of this type shall be said to be in quaternary canon- 
ical form. Just as a quaternion can be associated with a 4-tuple of reals as 



al + bi + cj + dk <=> (a,b,c,d), likewise the canonical linear function can be 
associated with a 4-tuple of quaternions as 

Aq + Bqi + Cqj + Dqk <^> {A,B,C,D} . 

The four-tuple of quaternions will be used as a shorthand notation for the 
canonical form. 

The linear sum of two such functions, say fi and f 2 , is immediately reducible 
to canonical form, i.e., 

fa (?) = fl (q) + h (q) = A 3 q + B 3 qi + C 3 qj + D 3 q k 

where 

A a =A 1 +A 2 , B 3 =B l +B 2 , C 3 = d+C 2l D 3 = D X +D 2 . 

However, functional composition of two such functions naturally expand into 
sixteen terms, which require reduction back to four terms. The number of terms 
in multiple compositions expand exponentially. The second reduction method 
provides a simple composition rule which automatically yields a quaternary 
form, provided the two functions being composed are already in this form. 

The next section provides preliminary material to introduce notation and 
nomenclature used in this work, as well as, matrix and involution operations 
useful in deriving the main results. Section[3]presents the method of reduction to 
quaternary form using either equivalent matrix operations or simpler involution 
operators. The matrix-based version provides a numerically efficient algorithm, 
while the involution-based version is better suited to symbolic analysis. Section 

01 presents the composition algorithm. 

Motivation behind this work is in the use of arrays of linear quaternion 
functions as filtering operations on quaternion valued data arrays. This allows 
for complex geometric operations on three- and four-dimensional data sets. Only 
after the basic manipulation of systems of such functions are detailed can the 
more difficult problem of linear quaternion filters be addressed. 

2 Preliminaries 

2.1 Quaternions 

The quaternions, H, are a four dimensional, non-commutative division algebra 
over the reals, R. Multiplication is defined by the hypercomplex operator rules 

i 2 = j 2 = k 2 — ijk = — 1. 

The typical quaternion is written in Cartesian form as 

P = Vo + P\i + P23 + P^,k 

where Po,i,2,3 S R- 

The quaternion conjugate, p, is defined as negating the non-real part, i.e., 
P = Po - P\i - P23 - Pzk 
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2.2 Quaternion-Matrix Homomorphism 



The equivalence between quaternion and 4x4 real matrix multiplication was 
first established by Peirce [3J, later verified by Spottiswoode g], made explicit 
by Cayley [51 [7] and later in Tait's book [5]. Walker 8 provides a concise 
introduction and Ickes [9] uses this equivalence to provide fast numerical com- 
putations on linear quaternion equations. 

These quaternion matrix representations are not unique. Farebrother et al. 
[TO] identified 48 different representations, not considering all automorphisms 
of the quaternion algebra. Later, Ducati |llj divided these representations into 
two groups based on the left and right action of the hypercomplex operators. 
Mackey [12] used this concept to solve the eigenvalue problem. Fassbender et 
al. [13] used them to solve the Hamiltonian eigenproblem. DeLeo and Rotelli 
[T"4"] used the same concept applied to complex 2x2 matrices. 

Using the hypercomplex operator product rules, the quaternion product of 
two arbitrary quaternions, q,p G H, written in Cartesian form as 

P = Po + Pii + P2.j + Psk 
q = q + qii + q 2 j + 93& 

is 

?~0 = P090 - Pl9l - P292 - P393 
ri = Piqo + Poqi - P3<72 + P293 
T2 = P2qo + P092 + P3qi - Pl93 

^3 = P3<7o + Po<?3 - P2qi + P\q2 
where r = pq = ro + r\i + rij + r^k. Gathering terms into matrix- vector 
notation we obtain two distinct matrix-vector products: the standard form 



?"0 




Po 


-Pi 


-P2 


~P3 




qo 


T\ 




Pi 


Po 


~P3 


P2 




qi 


T2 




P2 


P3 


Po 


~Pl 




92 


r 3 . 




. P3 


-Pi 


Pi 


Po 




. ?3 . 



(3) 



which preserves the order of the product, and the transmuted form 



ro 




9o 


-9i 


-92 


-93 




Po 


ri 




9i 


9o 


93 


-92 




Pi 


T2 




92 


-93 


9o 


9i 




P2 


r 3 . 




. 93 


92 


-9i 


9o 




. P3 . 



(4) 



which reverses the order. 

Following Ickes [SJ, let the standard matrix form be denoted [p] and the 
transmuted matrix form as [q] , where the 'f signifies transposing the lower- 
right 3x3 sub-matrix. 

In the quaternion product r — pq, if we regard p as a left-handed operator 
on the state variable q, then the p-operator is encoded as a 4x4 real matrix and 
the g-state as a real 4x1 column vector. Conversely, if q is regarded as a right- 
handed operator on the state variable p, then the operator is again encoded as 
a real matrix (transmuted to the left) and the state as a real vector. 
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2.3 Quaternion Operators 

Since left and right quaternion multiplications are unique, we can view various 
combinations as independent operators. Following De Leo and Rotelli [2], one 
can define a &ar-operator notation to provide a shorthand representation for 
these combinations. The &ar-operator, denoted ei|e2, is defined as: 

ei|e 2 (q) -> e 1 (q) e 2 , where e n e fc} 

so there are sixteen such operators as listed in tabular form below: 

1|1, l\i, l\j, l\k, 

i\l, i\k, 
il 1 , 

k\l, k\i, k\j, k\k. 

These quaternion operators are consistent with the standard hypercomplex 
operators and have similar operator formulae 

(i\l) 2 = (i|l) 2 = (k\lf = (i\l) (j|l) (k\l) = -l|l 
(l|i) 2 = (l|i) 2 = (l|fc) 2 = (l|i)(l|j)(l|fc) = -l|l 

so that various operator equations can be manipulated directly, e.g., (i\j) (j\k) = 
(ij\kj) = - (k\i). 

Using the standard and transmuted matrix homomorphisms given in the 
previous section, i.e., ei|e2 — > [ei] , each can be written as a matrix. For 
example, 



m _ [ k ] [iy = 









+1 














-1 


M 














-1 









The basic quaternion operator formulae can be extended by forming weighted 
linear combinations. Linearity of the quaternion operations implies that for four 
arbitrary quaternions, A, B,C, D 6 H, written in component form as 

A = a Q + a>ii + a 2 j + a^k, 
B = fa + Pi* + Paj + I3 3 k, 
C = 70 + 71 i + 72J + 73^, 
D = So + 5\i + S 2 j + Ssk, 

when combined with four operators 1|1, l|j, and l|fc, respectively, form the 
four extended operators 

(A\l) = a (1|1) + ai (i|l) + a 2 (j|l) + a 3 (k\l) , 
(B\i) = + t + (3 2 + p 3 (k\i) , 

(C|i) = 70 (i|i) + 7i (*|i) + 72 (ili) + 73 (k\j) , 

(D\k) = S (l\k) + 8 1 (i\k) + S 2 (j\k) + S 3 (k\k) . 
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In the standard case the quaternions ±1, ±i, ±j, and ±k form a non- 
abclian group of order eight, with multiplication as the group operator. In 
the quaternion-operator case, these become ±1|1, ±l|i, ±1|&, ±i|l, ±j|l, 

and ±fe|l giving a 14th order group with composition being the group operator. 

2.4 Quaternion Involutions 

A function, /, is an involution^ if it has the following properties [16] : 

1. is it's own inverse: /(/(?)) = q, 

2. is linear: f(q 1 + q 2 ) = f(q{) + f(q 2 ), and 

3. is anti-homomorphic: / (qiq 2 ) = /(<72)/(<7i)- 

A function is an anti-involutior^ if the third property, above, is homomor- 
phic, i.e., f(qiq 2 ) = /(tfO/faa)- 

Using the quaternion bar-operators in conjunction with conjugation provides 
a generalized conjugate, 

f -> -ege, ee{i,j,/c} 

for any gGi, which are involutions under these definitions |15J . The conjugate 
of these functions, 

F^-ege, e£{i,j,£;} 

are anti-involutions [17j . 

For an arbitrary quaternion q = qo + iqi + jq 2 + kq% the following identities 
hold in terms of involutions 

q\ = \ {q l + q) % 

q2 = \(q 3 +q)] (5) 

93 = 2 (q k + q)k 

which are analogous to the formula for the scalar part of a quaternion using the 
standard conjugate 

qo = | (q + q) (6) 

except they extract the remaining real-valued component-parts of any quater- 
nion. 

Finally, it should be pointed out that all of these functions (generalized 
conjugates, and component-part functions) are readily placed in quaternary 
form. This is achieved by first expressing the standard quaternion conjugate in 
involutory form, i.e., 

q = \ (-q + iqi + iqi + kqk) . 

Hence its canonical four-tuple is q = { — h, \, | , \ }. Proof of this involutory 
form is given in [TT] . 

1 This is also known as an involutory anti- automorphism. 
2 Likewise, this is also known as an involutory automorphism. 
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3 Reduction to Canonical Form 



3.1 Matrix Method 

At its core, the matrix method is the process for taking an arbitrary operator 
matrix and converting it to the sum of four operator matrices, each matrix 
equivalent to an operator quaternion. 

The matrix-quaternion homomorphism of section 12.21 suggest that the indi- 
vidual terms in equation |T]) can be converted into matrix form as 




where q denotes the state- vector form of quaternion q. 

If a numerical reduction to the problem is desired, then the following matrix 
equation could be used 

/(?)- |ek] 

Even though the resulting 4- vector represents a quaternion, the intermediate 
matrices [m p ] [n p ] cannot be decoded into a equivalent single quaternion since 
in general it does not have the correct pattern of elements. A general 4x4 
real matrix has sixteen degrees of freedom, whereas a quaternion has only four 
(whether or not it is in matrix form). So arbitrary real matrices cannot be 
written as a single quaternion. However, it can be written as the direct sum of 
four matrices, each matrix equivalent to a quaternion. 

Let A, B,C,D e H, be defined as 

A = a + a\i + a 2 j + a 3 k, 
B = [3 o + 0ii + 023 + fok, 
C = 70 + 7i i + 12] + Jak, 
D = S Q + 5ii + S 2 j + S 3 k, 

then each of the above quaternions contains four degrees of freedom. Hence the 
combined equation 

(A|l) + (B\i) + (C\j) + (D\k) 

has a total of sixteen; the number necessary to encode an arbitrary 4x4 matrix. 
Writing each quaternion operator in matrix form as 

-OL\ —Ol2 —C(3 

oto —as a 2 
as a —ot\ ' 
-ol-x a\ ao 



(B\i) -> [B] [i] T = 



-Pi 


-0o 


03 


-02 


A) 


-01 


-02 


-03 




-02 


01 


0o 


-02 


-03 


-0o 


01 
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(A\l) - [A] [l] t = 



Oil 

a 2 
as 



-72 


-73 


-70 


71 


-73 


72 


-7i 


-70 


7o 


-7i 


-72 


-73 


7i 


7o 


-73 


72 


-s 3 


s 2 


-Si 


So 


s 2 


S3 


So 


-Si 


-Si 


-So 


S3 


-s 2 


So 


-Si 


-s 2 


-S3 



(C\j) - [C] [jf = 



and 



(D\k) - [D] 



then the combined operator in matrix form is 
(A\l) + (B\i) + (C\j) + (D\k)^ 

®0 ~ Pi ~ 72 - 63 -«1 - Po - 73 + ^2 

ai + Po ~ 73 + ^2 ao - /?i + 72 + <5 3 
"2 + P3 + 7o - <5i a?3 - P2 - 7i - <5o 

"3 - /?2 + 7l + <5o -Q!2 - P3 + 10 - Si 

— «2 + /?3 — 70 — #1 -CK3 - P2 + 71 — <5o 

-a 3 - /?2 - 7i + <5o a 2 - /3 3 - 70 - <5i 
ao + /?i - 72 + S3 -ai + Po - 73 - S 2 
oti - Po - 73 - ^2 a + /3i + 72 - S3 

Setting this equal to an arbitrary matrix, R, with elements {r i3 }, and solving 
for the components of the quaternions one finds 



and 



D 



C = 



D = 



ao 




+r-ii + r 22 + r 33 + r 44 


U\ 


1 


+^21 - ri2 + r-43 - r 34 


a 2 


~ 4 


+^31 - ^"13 - ^42 + r-24 






_ +^4i - n 4 + r 32 - r 23 


" Po ' 




+r 2 i - r i2 - r 43 - r 34 


Pi 


1 


-ru - r 22 + r 33 + r 44 


P2 


_ 4 


-r 4 i - ri 4 - r 32 - r 2 3 


. P* . 




+r\3i + ri3 - r 42 - r 2i 


70 




+r 3 i - ri3 + r i2 - r 2A 


7l 


1 


+ru + ri4 - r 32 - r 23 


72 


~ 4 


-ru + r 2 2 - r 33 + r 44 


_ 73 




-r-2i - r-i2 - r-43 - r- 3 4 


" So ' 




+r 4 i - r-14 - r 32 + r 23 


Si 


1 


-r-31 - r-13 - r- 42 - r 2i 


s 2 


~ 4 


+r-2i + r-12 - r-4 3 - r-34 


S3 




-r-n + r-22 + r 33 - r 44 



(7) 



(8) 



(9) 



(10) 



This is not as difficult as it appears since the problem decomposes into 
four sets of four unknowns with each set independent of the others, e.g., the 



7 



coefficient set {ao, 72, £3} is solved using only the diagonal terms in both 
matrices. 

The steps in reducing equation |T]) to quaternary form are: 

1. Map coefficients m p and n p to matrix-space using equations ([3]|4]) 

trip -> [mp] , n p -> [n p ] f . 

2. Construct single matrix R as 



R = ^[m p ] [rip] 1 . 

p=l 

3. Construct equivalent quaternion factors using equations (T7l ll0[) 
fl-(A|l) + (B|i) + (C|j) + (D|fc). 

3.2 Involution Method 

The matrix equivalent method outlined above is useful for numerical computa- 
tions, however, it requires manipulation at the component level of each quater- 
nion, m p and n p . The involution method, outlined in this section, provides 
an alternative solution in terms of m p and n p directly, without these explicit 
references, making it more suitable for symbolic analysis. 

Starting with equation |T|). let each n p be expanded into component form as 

n p = w p + ix p + jy p + kz p 

then the full equation can be expanded as 

p 



f (?) = ^2 m p q (w p + ixp + jy p + kz p ) 
P =i 

p p 
= m p w p q + m p x p q i 
p—i p— 1 

p P 

+ ^2 mpVpQ 3 + ^2 m p z p1 k 



p= 


1 






' p 




" p 






q + 


^ TTl p X p 


qi 


J3=l 




.P=l 




■ p 




' P 





.p=i 



. P =i 



q k. 
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Now let 

p 

A = rripWp, B 
P =i 
p 

C = m p y p , D 
p=i 

so that 

f (q) = Aq + Bqi + Cqj + Dq k. 

The problem with this formula is that we are required to decompose quaternions 
n p into component form, whereas we are looking for a formula that does not 
require this step. The answer to this conundrum is to apply the involution 
formulas ([5][B]). Namely, let 

w p = \ {n p + n p ), x P = \ {np + n p ) i, 

Dp = \ (n p + fip) j, z p — \ {™> P + %>) ^> 

so that 

p 

A=\^2m p (np+fip), 
p=i 
p 

B = | m p (n p + n p ) i, 
P =i 
p 

p=l 
p 

D = i 2_. m p { n p + n p) k- 

p=l 

4 System Reduction 

One could use the methods of the previous section to reduce compositions of two 
linear functions, i.e., where f 2 (q) o f x (q) = f 2 (/i (q)). What is not apparent 
is that functional composition can be reduced in a more elegant way. This 
alternative lends itself to better insight into what happens in these composition 
operations. 

Using bar-operator notation, composition is reduced to ordered multiplica- 
tion, hence one can build a composition table showing how components of the 
composed functions are reduced. The composition table for f 2 ° /i is given below 



/ 2 °/i 


Ai\l 


B x \i 


Ci\j 


Di\k 


A 2 \l 


+A 2 A ± \l 


+A 2 B x \i 


+A 2 C 1 \j 


+A 2 D x \k 


B 2 \i 


+B 2 Ai\i 


-B 2 B X \1 


-B 2 d\k 


+B 2 Dt\j 


c 2 \j 


+C 2 At\j 


+C 2j Bl|fc 


-C 2 Ci\\ 


-C 2 Dx\i 


D 2 \k 


+D 2 A ± \k 


-D 2 B x \j 


+D 2 d\i 


-D 2 D X \\ 



E 

p=i 

p 



p=l 
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Collecting like terms, we find that 



h h (q) = A 3Q + B 3 qi + C 3 qj + D 3 q k 



where 



^3 

B 3 
C 3 
D 3 



A 2 A X - B 2 B 1 - C 2 Ci - D 2 D U 
A 2 B 1 + B 2 Ai - C 2 Di + D 2 d, 
A 2 C X + B 2 D 1 + C 2 A X - D 2 B X , 
A 2 D X - B 2 d + C 2 B X + D 2 Al 



Careful examination of the 4-tuple resulting from a product of two quater- 
nions with the 4-tuple resulting from the functional composition of two canoni- 
cal equations reveals that they are transmuted versions of the same operational 
formula. This is a consequence of the bar-operators having a similar operator 
formula as the standard hypercomplex operator formula. If the canonical form 
had been defined with the full quaternions on the right, e.g., 



then the composition formula would have been identical to the quaternion mul- 
tiplication formula. 

5 Conclusions 

The goal of this work is the reduction and manipulation of systems of linear 
quaternion functions. To this end two results have been provided. First, even 
though a general linear quaternion function may contain an arbitrary number of 
terms, it can always be reduced to at most four terms. This reduced function, 
called the quaternary canonical form of the function, is completely specified with 
a four-tuple of quaternions. Second, series (cascaded composition) and parallel 
(weighted-sum) combinations of such canonical functions can be written directly 
in canonical form without recourse to the previous reduction method. 
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